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Abstract 

In this letter, we introduce a geometric model to explain the origin of the observed shallow levels in 
semiconductors threaded by a dislocation density. We show that a uniform distribution of screw dislocations 
acts as an effective uniform magnetic field which yields bound states for a spin-half quantum particle, even 
in the presence of a repulsive Coulomb-like potential. This introduces energy levels within the band gap, 
increasing the carrier concentration in the region threaded by the dislocation density and adding additional 
■ recombination paths other than the near band-edge recombination. 
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I. INTRODUCTION 



Dislocations in semiconductors are a nuisance for device developers since they introduce states 
which trap electrical charge, reducing the number of available carriers. By getting charged they 
introduce electric fields, locally affecting device performance by scattering the electrons. Their 
states act as non-radiative recombination centers for electrons and holes, therefore reducing the 
efficiency of opto-electronic devices. They may also cause current leakages by jumping of the 
electrons from state to state. In this work we introduce a geometric model to describe the effect 
of a density of screw dislocations on a spin-half charged particle. We show that torsion, associated 
to the defect distribution, acts as a magnetic field on the particles, giving rise to bound Landau- 
like levels even in the presence of a repulsive electric field. Since the magnetic- like effects due to 
torsion are insensitive to the signal of the electric charge, depending on the location of the Fermi 
level one might have either bound holes or electrons in these states. This way, a suitable applied 
magnetic field may be used to partially cancel the effects of torsion freeing some of the trapped 
charge carriers from the bound states. 

The concern with the influence of dislocations on charge carrier mobility in semiconductors is 
not new Shockley [2J was the first to suggest that dangling bonds in the dislocation core may 
act as traps (deep levels). Since these earlier studies a large number of articles have been published 
on the effects of dislocations in the properties of semiconductors. For reviews, see jjj. Besides the 
dangling bond traps, electrons and holes may be trapped in more extended states (shallow levels) 
due to the elastic deformation field of a dislocation or of a density of dislocations. This is the 
subject of this article, where the deformation field enters as a geometrical object, torsion. 

In the elastic continuum, a torsion field appears linked to the strain and stress introduced by 
a topological defect distribution 5, (J. Katanaev and Volovich jf| have shown the equivalence be- 
tween the continuum theory of defects in elastic solids and three-dimensional gravity with torsion 
(that is 3D Riemann-Cartan geometry). The deformation introduced by the defect is described 
geometrically by a metric which corresponds to a particular solution of the 3D Einstein-Cartan 
equation [J, 6, 7]. Studies of the influence of topological defects on quantum systems using the 



Katanaev- Volovich approach have been made in quantum scattering 
ativistic scalar particle [9( , the self- 



tral particle 
quantum ring 
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Berry's phase 



8( , Landau levels for a nonrel- 
adjoint extension method 1.0 . Landau levels for a spin-half neu- 



121 ]. holonomic quantum computation [13J, and a two-dimensional 



141 ] . Other studies of topological defects related to a torsion field, from the classical 



point of view, are geodesies around a dislocation [151 ] and torsion effects on electromagnetic fields 
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In this letter, with the purpose of better understanding the shallow levels that appear in semi- 
conductors threaded by screw dislocation densities, we study the behavior of a spin-half charged 
quantum particle under the influence of a Coulomb-like potential in an elastic medium containing 
a uniform distribution of screw dislocations. We show that the uniform distribution of screw dislo- 
cations plays the role of an effective uniform magnetic field, and yields bound states for a spin-half 
quantum particle even when the Coulomb-like potential is repulsive. Our results suggest that one 
might use a magnetic field to destroy the dislocation density bound states and therefore minimize 
its effects on the charge carriers. 



II. GEOMETRIC MODEL 

n 

As shown in [17j] a spinless quantum particle moving through a uniform distribution of screw 
dislocations attains bound states analogous to the Landau levels of a charged particle in a uniform 
magnetic field. Even though the torsion associated to the defect distribution does not affect charge 
it does couple to spin. Electrons and holes are therefore undistinguishable from the point of view 
of torsion. In order to get as close as possible to the problem of carrier motion in the presence of 
a distribution of dislocations, we extend [17!] to include spin and a Coulomb-like potential. The 
geometry corresponding to a uniform distribution of parallel screw dislocations is described by the 
line element 17J in cylindrical coordinates 

ds 2 = dp 2 + p 2 dip 2 + (dz + Q,p 2 dip) 2 , (1) 

where = b z ^, with A being the area density of dislocations and b z the Burgers vector. The 
z-axis was chosen to lie parallel to the screw dislocations. 
With the choice of the nonholonomic 1-form basis 

1 = dp; 6 2 = pdtp; 6 3 = dz + Qp 2 dip, (2) 

the line element ([1]) takes the form ds 2 = 9 a 9 b 5 a b. The non-coordinate basis ([2]) is related to the 
holonomic (flat space) basis 

dx 1 = dp; dx 2 = pdip; dx 3 = dz, (3) 

by the transformation 6 a = e a ^(x)dx^. From now on we reserve latin indices to objects writ- 
ten in terms of the nonholonomic basis and greek indices for the ones written in terms of 
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the holonomic basis. The transformation matrix e a ^ (x) is known as triad. It follows that 
ds 2 = g^ v (x) dx^dx v = e a ' (x) dx fJ- e b u (x) dx u 5 a i, and therefore (x) = e a (x) e b v (x) 5 a t,. That 
is, the geometry is encoded in e a ^ (x) e b u (x) . Furthermore, the triad has an inverse given by 
dx^ = e^ a (x) 8 a , where both, triad and its inverse, satisfy the relations: e" (x) e^ b (x) = 5 a b and 
e\ (x) e\ (x) = n. 

n 

Following the formulation of the spinor theory in curved spacetime [181 ]. we have that the par- 
tial derivative becomes the covariant derivative of a spinor whose components are defined as 
V M = <9 M + f M (x), where f M (x) = | oj^ab {%) ^ ab is called the spinorial connection, uj^ a b (%) is a 
connection 1-form related to the curvature of the manifold, and S a6 is defined for two-spinors 
by S ab = | [cr a ,cr b ]. We denote cr° = / as the 2x2 identity matrix, and the matrices a 1 
correspond to the usual Pauli matrices. In the presence of torsion, the expression of the co- 
variant derivative of a spinor changes [ijj]. In this case, the spinorial connection is defined by: 
r /t (x) = | [oo^ab i x ) + Kfiab (%)] ^ ab ■ The connection 1-form (x) is related to the contortion 
tensor via the expression [19]: 

(x) = K/3 Ufl (x) e u a (x) e P b (x) - e u b (x) e p a (x) . (4) 

Moreover, the contortion tensor Kp u/1 (x) is related to the torsion tensor (x) via 







T f/i (s) T v ^ (x) u (x) 



Note that the torsion tensor is antisymmetric in the last two 
indices, while the contortion tensor is antisymmetric in the first two indices. Following these 
definitions, we can also represent the torsion tensor in terms of three irreducible components: the 
trace vector T M = ^ the axial vector S a = e a/3u ^ T^^^ and the tensor qp UIJL which satisfies 
the conditions: q a = and qp Ufl = 0. In this way, the torsion tensor can be written as 

Tpufi = I (Tu9pn ~ T ^9/3u) - I epviey + qpufi, and the connection 1-form flU) can be defined in 
terms of these irreducible components [19 [.By writing the torsion tensor in terms of the irreducible 



components, it has been shown in Ref. 



191 that the axial 4-vector S 1 ^ couples to spinors. Hence, 

n 



the Schrodinger-Pauli equation in the presence of curvature and torsion is given by [201 ] 

.dtp 1 



dt 2m 



iP + \o-S^ + V{p) if>, (5) 



where the vector E is defined in such a way that its components are given in the local reference 
frame by Ek = \ a 3 e' f k (x) — | S° o^. As discussed in [jj]], the components of the vector a can be 
considered as internal degrees of freedom, that is, \a corresponds to the spin of the particle. In 
this way, the coupling between the 4-vector and spinors gives rise to the term | a ■ S in Eq. 
©, which is called the spin-torsion coupling and is analogous to the Zeeman spin- magnetic field 



coupling which introduces a splitting of each of the states into a pair, one for spin up particles 
and the other for spin down particles. Also, we have a term coupling the linear momentum to the 
spin, — jSq<j ■ p, analogous to helicity. In order to solve the Schrodinger-Pauli equation ([5]), we 
need to note that both connections 1-form oj^ b (x) and K a b (x) can be obtained by solving the 
Cartan structure equations 21] T a = dO a + uj a b f\9 b , where the operator d is the exterior derivative, 



the symbol A means the wedge product, uj a b = u' a b (x) dx^ is the spin connection 1-form, and 



T a _ 1 T a^ dx fM A dx u ig the torsion 2-form. 



Now, by solving the Cartan structure equations for the triads given in ([2]), we obtain u 2 ( x 



2 

<P 1 



-1 and T 3 = 2Qpdp A dip. Hence, we obtain [22j only one non-null component of the 



axial 4-vector S^, which is S° = —40. Furthermore, let us consider a Coulomb-like potential given 
by: 

V(p) = f - = ±\A i (6) 

where / is a constant. Note that the plus (minus) sign in Q means that the Coulomb-like potential 
is repulsive (attractive). We shall see that, bound states can be achieved for either sign of the 
Coulomb-like potential ([6]) due to the influence of the uniform distribution of screw dislocations. 
This way, by considering the spin being aligned with the symmetry axis of the screw dislocations 
(z-axis), the Schrodinger-Pauli equation ([5]) becomes 
.dip 1 



dt 2m 



d 2 ld_ 1 d 2 ^ d 2 . 2 2 , d 2 
dp 2 p dp p 2 dp 2 dzdp ^ dz 2 



2m p 2 dip 

(7) 

ia 3 d^ i . n 2 ( / , 

2m dz 8mp 2 8m p 



We can see in Eq. (|7|) that ip is an eigenfunction of cr 3 , whose eigenvalues are s = ±1 and the 
Hamiltonian of Eq. ([7|) commutes with the operators [23] J 2 = — id v and p z = —id z , thus, we can 
write the solution of equation (|7|) in terms of the eigenfunctions of the operators J z and p z , that 
is, ip s = e~ l£t e l ^ l+ ^ ip e lkz R s (p), where I = 0, ±1, ±2, . . . and A; is a constant which corresponds to 
the momentum in the z-direction. We take k > since, as we will see below in equation (jlOp . for 
k < the minus sign of the exponent of the Gaussian function becomes positive and we no longer 
have bound states. This asymmetry is due to the choice of the Burgers vector orientation. In this 
way, substituting this general solution into the second order differential equation (J7J), we have 



where we have defined 7 S = I + A (1 — s). By making a change of variables given by £ = \/0A; p, 
we can rewrite the radial equation ([8]) in the form: 



d 2 R s 1 di? s 7s 2 



-Rs - C 2 Rs 7^f~r + Rs - °> 



(9) 



dc 2 c dc e 

where /3 S = 2m £ — (k + s §) 2 — — 7 S and /' = 2m/. The solution of © should be given in 
such a way that the wave function is regular at the origin. In this way, we can write the solution 
of (l9l) in the form: 



R s (C) = e-^(^H S (() 



(10) 



Hence, substituting (|10|) into (|9~]), we obtain the following second order differential equation: 



+ 



+ 2C 



+ 



/' 



H s = 0, 



(11) 



where we have defined in (jlip the parameters: a = 2 |7 S | + 1 and g = ^ — 2\j s \ — 2. The function 
which is solution of the second order differential equation (jlip is the Heun biconfluent function 



24. 



25fl: 



H a (0 = # 



91 I n & ^ 



OA;' x/oT 



(12) 



In order to look for bound states, for both signs of the Coulomb-like potential given in ([B]), we 
use Frobenius method [3] to write the solution of Eq. (fTTI) as a power series expansion around 



the origin, that is, 



m=0 



Substituting the series (113f) into (jlip . we obtain the recurrence relation: 

/' a m+ i (5 - 2m) 

«m+2 



(13) 



(14) 



/fifc (m + 2) (m + a + 1) (m + 2)(m + a + l) 
By starting with ao = 1 and using the relation (I14p . we can calculate the first three coefficients 
of the power series expansion (113p . Then, we have 



01 



«2 



03 



/' 2 



3 



2a (a + 1) OA; 2 (a + 1) ' 
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(15) 



1 



^7 6a (a + 1) (a + 2) 



OA: J 6a (a + 1) (a + 2) 



(g 



QkJ 3a (a + 2) 



Note that the recurrence relation (|14p is valid for both signs of the Coulomb-like potential, that 
is, we can specify the signs of the Coulomb- like potential by making / — > ± |/| in (|14p . Hence, in 
order to obtain finite solutions everywhere, which represent bound state solutions, we need that 
the power series expansion (|13p or the Heun biconfluent series become a polynomial of degree re. 
Through expression (|14p . we can see that the power series expansion (|13p becomes a polynomial 
of degree n if we impose the conditions: 

g = 2n and a n+ \ = 0, (16) 

where re = 1, 2, 3, . . .. Prom the condition g = 2n, we obtain the expression for the energy levels of 
the bound states. Thus, we have 



1 

2m 



k + s- 



V 2 



(17) 



where the angular frequency is given by w nj / jS = — . On the other hand, the condition a n+ \ = 
allows us to obtain a expression involving the angular frequency and the quantum numbers n, I and 
s. Observe that we have considered k being a positive constant, thus, we can choose any value of k 
in such a way that the condition a n +i = and write k = k n j s . We should note that writing k in 
terms of the quantum numbers re, I and s, that is, k = k n i s , it does not mean that k is quantized. 
Writing k = k n j s means that the choice of the values of k > depends on the quantum numbers 



re, I and s in order to satisfy the condition a n +i = 
as i,s = —^ LL± - In the following, we c 
and n = 2. In this way, for re = 1, we have 



32 ]. Then, we can write the angular frequency 



as uo n i s = *' 3 . In the following, we calculate the values of the angular frequency for n = 1 



^-^^"2^(217,1 + 1)' (18) 



and for re 



nk 2 ,i, s f' 2 nQ , 

W2 ' , ''- _ ^-4m(4| 78 |+3)- (19) 

Hence, we have shown that bound state solutions of the Schrodinger equation ([7]) can be achieved 

for both attractive and repulsive Coulomb-like potential ([6]) in the presence of a uniform distribution 

of screw dislocations. Due to the asymmetry introduced by the direction of the Burgers vector the 

bound states appear only for motion in the positive z-direction. The origin of the bound states 

is the coupling between the angular variable ip and the linear variable z introduced by the screw 

dislocation density. If there is no motion along the z-axis, there is no coupling and hence, no bound 

state. This can be seen by inspecting equation (|10p . With k = the Gaussian function, which 

assures we have bound states, becomes unity. 



III. CONCLUSIONS 



In this letter, we studied the behavior of a spin-half quantum particle under the influence of a 
Coulomb-like potential in an elastic medium containing a uniform distribution of screw dislocations. 
We showed that the uniform distribution of screw dislocations plays the role of an effective uniform 
magnetic field, and yields bound states for a spin-half quantum particle even when the Coulomb- 
like potential is repulsive. We obtained a finite solution of the radial equation by imposing that the 
Heun biconfluent series becomes a polynomial of degree n, which is possible if the two conditions 
given by (fT6|) are satisfied. 

We saw that a density of screw dislocations may bind charge carriers in Landau-like levels in the 
transversal plane, while leaving them free to propagate along the direction of the Burgers vector, 
but not in the contrary direction. Application of a uniform magnetic field may selectively cancel 
this effect since the magnetic field couples to electric charge but the screw dislocation density does 
not. In other words, torsion, due to the dislocation density, puts both electron and holes to move 
in cyclotron orbits in the same sense. The magnetic field does the same with electrons and holes 
moving in reversed senses. Depending on the direction of the applied field either electrons or holes 
are freed from their bound states. 

Even though dislocations are generally seen as problems, since 1979 there have appeared pro- 
posals 



281430] of using them as active one-dimensional conductive channels for applications. There 
is even the possibility of having one-dimensional topologically protected conducting modes associ- 
ated to dislocations in topological insulators [3l| . We hope this work contributes towards uses of 
these "undesired" defects in semiconductor devices. Quoting George Bernard Shaw, "If you can't 
get rid of the skeleton in your closet, you'd best teach it to dance." 
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